NREP, transcriptionally upregulated by HIF-1α, aggravates breast cancer cell growth and metastasis by promoting glycolysis

Breast cancer (BC) poses a great threat to women’s health. Neuronal regeneration related protein (NREP) is a multifunctional protein that is involved in embryonic development, regeneration, and human disease. However, the biological function of NREP in tumors is rarely reported and its role in BC remains unknown. Bioinformatics analysis showed that NREP is highly expressed and closely correlated with poor survival in BC patients. Under hypoxic conditions, NREP was upregulated in BC cells, and this promotion was reversed by hypoxia-inducible factor HIF-1α suppression. Luciferase reporter system and chromatin immunoprecipitation assays confirmed that HIF-1α directly binds to the promoter of NREP to increase the transcriptional activity of NREP. NREP suppression inhibited cell proliferation, arrested the cell cycle at the G1/S phase, and promoted apoptosis and caspase-3 activity in BC cells. Suppression of NREP decreased the tube formation ability of HUVECs. In addition, NREP downregulation showed an inhibition effect on cell migration, invasion, and EMT of BC cells. In NREP overexpressed cells, all these changes were reversed. In vivo, animal experiments also confirmed that NREP promotes BC tumor growth and metastasis. In addition, NREP promoted cellular glycolysis and enhanced the levels of glucose consumption, ATP, lactate production, and glucose transporters expression in NREP-overexpressed BC cells. In summary, our results demonstrated that NREP could be transcriptional activated by HIF-1α, which may aggravate BC tumor growth and metastasis by promoting cellular glycolysis. This result suggested that NREP may play an essential part in BC progression.


INTRODUCTION
Breast cancer (BC) is a complex disease that displays both interand intra-tumoral heterogeneity [1] and is the most frequently diagnosed cancer and the leading cause of cancer-related death among women worldwide [2].It also has a complicated pathogenesis, including hormone disorders, genetics, environment, etc. Surgery, chemotherapy, hormone therapy, and radiotherapy are common treatments for BC clinical patients [3].However, metastasis and recurrence remain the main causes of poor survival in BC patients, even though significant advances have been made in the treatment of primary BC [4][5][6].Therefore, it is of great significance to explore the key factors during tumor progression, which may alleviate the clinical outcome of BC.
Neuronal regeneration related protein (NREP) is a conserved RNA-binding intracellular protein containing only 68 amino acids [7].NREP was first identified by Studler and colleagues [8] and is relatively conserved in many species [9], expressed in the brain, smooth muscles, regenerated tissue, and malignant glioblastomas [10].NREP plays an essential role in promoting wound closure [11].Reported studies showed that NREP promotes the proliferation and migration of fibroblasts [12][13][14], and the migration of hepatic astrocytes [15].Nevertheless, studies on the exploration of NREP in tumorigenesis is largely limited.In gastric cancer [10,16], NREP promotes cell proliferation, angiogenesis, and EMT process.In human glioma cells, NREP inhibition decreased cell migration [17].In lung cancer cell A549, NREP could promote cell cycle progression [18].These findings indicate that NREP may play an important role in tumor cell proliferation and migration.However, its role in BC is still unknown.
Tumor cells mainly depend on glycolysis to produce energy [19][20][21].Therefore, a targeted regulation of glycolysis may provide a therapeutic direction for BC treatment.Reported studies displayed that NREP could activate TGF-β1 [22], which promotes the malignant progression of tumor cells by inducing glycolysis through metabolic reprogramming [23,24].Also, NREP activates Rho GTPases family proteins like RhoA, Rac1, and RalA [25][26][27].It has been proved that inhibition of Rac1 inhibits glycolysis and proliferation of BC cells [28].Interestingly, BC is proved to have a considerably lower oxygen partial pressure than that in normal breast tissue [29], implying a hypoxic microenvironment of BC.Under hypoxia conditions, hypoxia-inducible factor HIF-1α promotes glycolysis of tumor cells, including BC.In addition, HIFs could activate the transcription of downstream factors, thereby participating in tumor progression [30].Importantly, bioinformatics analysis showed that NREP is upregulated in hypoxia BC cells than in normoxia.Therefore, we speculated that NREP may be upregulated by HIF-1α, and involved in cellular glycolysis progression of BC.
In this study, we aimed to investigate the biological function role of NREP in proliferation and metastasis of BC, and the potential up-/down-stream regulatory mechanisms involved.

NREP is upregulated and predicts poor prognosis value in BC clinical samples
To explore the genes that may play an important role in promoting BC progression, we used 9 public available BC databases from GEO to screen the differentially expressed genes (DEGs).Venn diagram showed that there were 22 overlapping upregulated DEGs (Fig. 1A).The specific expression level of these 22 genes in each database was displayed as the heat map (Fig. 1B).Then, we used functional enrichment analyses to explore the function of these genes.For KEGG (Fig. 1C), these genes were related to the terms named "chemokine signaling pathway", and "pathways in cancer"; and for GO (Fig. 1D), these genes were related to the terms named "positive regulation of cell migration", "apoptotic process", and "identical protein binding".Among these genes, we noticed NREP, which exhibited a relatively higher expression level in all of these 9 databases.Compared with the normal tissues, tumor tissues exhibited a significantly higher expression level of NREP in these 9 databases (Fig. 1E).Besides that, NREP was also upregulated in BC tumor tissues as compared with the para-carcinoma and healthy tissues (Fig. 1F, analyzed by Breast Cancer Gene-Expression Miner v4.9).Survival curves showed that high-level NREP in BC patients displayed poor survival probability in overall survival (OS), metastasis-free survival (MFS), progression-free interval (PFI), and progression-free survival (PFS) (p < 0.05, Fig. 1G).Besides that, in BC patients, NREP exhibited higher expression in the tumor tissues when compared to the para-carcinoma tissues (Fig. 1H).Therefore, we suspected that the enhanced level of NREP may take a great part in BC tumor development.

NREP is transcriptionally upregulated by HIF-1α in BC cells
From GEO databases (GSE3188 and GSE111259), we found that NREP is upregulated in hypoxia BC tumor cells when compared to normoxia (Fig. 2A).Depletion of HIF-1α in BC cells exhibited an inhibition level of NREP (data from GSE3188, Fig. 2B).Then, we tested the expression of NREP in BC cell lines SK-BR-3, MCF-7, MDA-MB-468, and MDA-MB-231 (Fig. 2C).MDA-MB-468 and MDA-MB-231 showed a more obvious change between normoxia and hypoxia conditions.Thus, these two cell lines were chosen to achieve NREP up-or downregulation.Consist with the results from GEO databases, HIF-1α and NREP showed an increased mRNA (Fig. 2D, E) and protein (Fig. 2F) expression in hypoxia than normoxia, while downregulation of HIF-1α reversed this enhancement in BC cell lines.These results confirmed that hypoxia affects the expression of NREP in BC cells.
To verify if NREP could be directly regulated by HIF-1α, we used a dual-luciferase reporter system and ChIP-PCR assay for confirmation.In Fig. 2G, hypoxia caused a significantly increased transcriptional activity of the NREP promoter, while the downregulation of HIF-1α decreased this enhancement.In Fig. 2H, serial deletion results revealed the region between −735 to −383 bp (site 2) was essential to NREP-induced expression of luciferase reporter (significantly decreased luciferase activity).However, the regions between −1733 to −735 bp (site 1) and −383 to −12 bp (site 3) exhibited almost the same relative luciferase activity, implying that the putative HIF-1α binding sites in these two regions were unimportant.Besides that, the region between −12 to +193 bp (site 4) still has a high relative luciferase activity, suggesting that this site may also be important.Therefore, we designed primers for sites 2 and 4 (Fig. 2I) and performed a ChIP-PCR experiment.In Fig. 2J, protein bands appeared at sites 2 and 4 after immunoprecipitation.
IgG used as a negative control exhibited no bands.This result confirmed that HIF-1α could directly bind to the promoter of NREP, thus upregulating its expression.
NREP promotes BC cell growth and angiogenesis NREP was down-or up-regulated in MDA-MB-468 and MDA-MB-231 cells, and mRNA (Fig. 3A, B) and protein expression (Fig. 3C, D) were verified.Then, we explored the effect of NREP on BC cell proliferation and apoptosis.Results showed that both in MDA-MB-468 and MDA-MB-231 cells, NREP knockdown suppressed cell viability (Fig. 3E, F), blocked the G1 to S phase transition in the cell cycle (Fig. 3G, H), decreased colony formation rate (Fig. 3I, J) and EDU positive cells (Fig. 3K, L), while all these changes were reversed in NREP overexpressed cells.In addition, NREP silencing increased the cell apoptosis rate (Fig. 3M) and caspase-3 activity (Fig. 3N) in MDA-MB-468 or MDA-MB-231 cells.Moreover, we investigated the effect of NREP on angiogenesis.HUVECs were cultured in tumor-conditioned medium (TCM) from NREP up-or down-regulated MDA-MB-468 and MDA-MB-231 cells.NREP silencing inhibited tube-like structure formation in HUVECs, while NREP overexpression induced it (Fig. 3O, P).

NREP promotes BC cell migration, invasion and metastasis
After that, we explored the effect of NREP on BC metastasis in MDA-MB-231 cells.Images and quantification data from the Transwell assay demonstrated that knockdown of NREP inhibited the numbers of migrated (Fig. 4A, B) and invaded (Fig. 4C, D) cells while overexpressing NREP promoted it.In addition, MDA-MB-231 cells exhibited striking morphological modifications.Upregulation of NREP triggered the transition from a cobblestone-like morphology to an elongated shape associated with increased cell scattering (Fig. 4E).Also, the E-cadherin expression was upregulated, and Ncadherin, SLUG, and MMP9 protein expression were decreased in response to NREP inhibition (Fig. 4F).Meanwhile, the mRNA expression of E-cadherin was upregulated (Fig. 4G), and N-cadherin was downregulated (Fig. 4H) with NREP knockdown.In NREP upregulated MDA-MB-231 cells, all these trends were reversed.
Then, we analyzed the effect of NREP on BC metastasis in vivo.In Fig. 4I, the metastatic nodules in the lungs of mice injected with NREP-knockdown cells were less than that of mice injected with control cells, while more metastatic nodules were apparent in mice injected with NREP-overexpression cells than of vector infected cells.Similarly, the H&E staining images of lung tissues (Fig. 4J) and quantification data of metastatic nodules (Fig. 4K) also showed the same result.

NREP promotes BC tumor growth in vivo
Based on these results, we confirmed the promotion effect of NREP on BC progression in vivo.Mice were injected with NREP-sh or -ove infected MDA-MB-468 cells.The images and volume changes of tumors were shown in Fig. 5A, B. Similar to the in vitro result, the knockdown of NREP decreased tumor growth, while overexpression of NREP showed an opposite result that promotes tumor growth.The decrease of NREP (Fig. 5C) and Ki67 (Fig. 5D) mRNA expression, fewer NREP-positive cells (Fig. 5E) and Ki67-positive cells (Fig. 5F), and an increase of red fluorescence of TUNEL images (Fig. 5G) were observed in NREP silenced cells formed tumor tissues.Besides that, suppressed CD31 mRNA level (Fig. 5H), decreased microvessel density (Fig. 5I), increased E-cadherin mRNA level (Fig. 5J), and decreased N-cadherin mRNA level (Fig. 5K) were also observed in NREP silenced cell formed tumor tissues.In NREP upregulated cells formed tumor tissues, and all these trends were reversed.

NREP accelerates cellular glycolysis in BC cells
After that, we sought to determine the potential role of NREP in the glucose metabolism of BC.As shown in Fig. 6A-C, either in MDA-MB-468 or MDA-MB-231, NREP regulation displayed the same trend.NREP knockdown decreased cellular glucose consumption, ATP level, and lactate production, while upregulation of NREP led to increased levels in BC cells.Then, we detected the protein expression of glucose transporters in BC cells.As exhibited, NREP inhibition had a lower protein level of GLUT1, HK2, PKM2, and LDHA in MDA-MB-468 (Fig. 6D, E) or MDA-MB-231 (Fig. 6F, G) cells.Conversely, in NREP upregulated cells, all of these genes were upregulated.Therefore, these results indicated that NREP may promote the growth and metastasis of BC by regulating cellular glycolysis.

RNA-seq and functional enrichment analyses of NREP silenced BC cells
After that, RNA-seq was used in NREP-inhibited MDA-MB-468 cells to explore the downstream factors and pathways of NREP.PCA plot (Fig. 7A, B) and hierarchical clustering plot (Fig. 7C) were performed to confirm the consistency and variance of samples.Volcano map exhibited all DEGs, including 685 upregulated DEGs and 1898 downregulated DEGs (Fig. 7D).Then, we performed the GO and KEGG functional enrichment analyses in all DEGs.The top 20 terms of GO were mainly enriched in "protein binding", "cytosol" and "cytoplasm" (Fig. S1A).As for KEGG, a total of 51 pathways were enriched (p < 0.05).The pathways that related to "Cellular Processes" and "Metabolism" were exhibited in Fig. 7E, and the others were shown in Fig. S1B.Then, the DEGs that enriched in "cell cycle" and "glycolysis/gluconeogenesis" pathways were focused.We confirmed the effect of NREP regulation on CCNA2 and CCNE1 expression, both of them played an important   role in the cell cycle G1/S transition.Downregulation of NREP suppressed the mRNA expression of CCNA2 and CCNE1 in MDA-MB-468 cells, while overexpressing of NREP promoted their expression (Fig. 7F).Besides that, based on RNA-seq results, we showed the mRNA expression of glycolysis-related factors.Downregulation of NREP suppressed the mRNA expression of LDHA, GLUT1, and PKM2 in MDA-MB-468 cells (Fig. 7G).

DISCUSSION
As a malignant disease threatening women's health, BC has always been the focus of the majority of studies in cancer.In this work, we found that NREP is highly expressed in BC tumor tissues and cell lines.HIF-1α as an upstream factor transcriptionally activated NREP and upregulated its expression.Highly expressed NREP promoted cell proliferation, angiogenesis, migration, invasion, and EMT process of BC cells in vitro and in vivo.Moreover, NREP promoted glycolysis of BC cells.Accordingly, we demonstrated that aberrantly expressed NREP may play a promoting role in tumor growth and metastasis in BC by promoting cell proliferation and glycolysis.
Hypoxia is an essential factor in cancer development [31][32][33][34], and this effect was related to the activation of downstream genes that are subject to HIF-1α [35][36][37].In BC cells, HIF-1α plays a critical role in stimulating the metastasis [37] of primary tumor to distant organs, which is closely related to patient mortality.Knockdown of HIF-1α abrogates this effect of hypoxia and significantly impairs the tumor-initiating ability of BC cells [31].These findings indicated that hypoxia is essential for BC tumor development.Therefore, after confirming that NREP is highly expressed in different BC databases and that its expression may be induced by hypoxia, we first verified whether NREP functions as downstream of HIF-1α.Consistent with online analysis, our results showed that knockdown of HIF-1α decreased the level of NREP in BC cells.Dual luciferase reporter and ChIP-PCR assay confirmed that HIF-1α could bind to the promoter of NREP, thus regulating its expression.Therefore, we referred that the enhanced level of NREP may take a great part in BC tumor development, and this process may be related to HIF-1α.
After that, we explore the effect of NREP on the malignant phenotype of BC cells in vitro and in vivo.So far, the function of NREP in gastric cancer has been explored, and they found that NREP promotes the growth and metastasis of gastric cancer [10,16].Similar to the reported study [10,16], upregulation of NREP promotes cell proliferation, angiogenesis, migration, invasion, and EMT process in vitro and in vivo.These findings indicated that NREP played an oncogene role in BC tumor growth and metastasis.In lung cancer cell A549, NREP could promote the G1/S phase transition of the cell cycle [18].Therefore, we examined the effect of NREP on the BC cell cycle.Results showed that NREP downregulation blocked the G1 to S phase transition in the cell cycle, while overexpressed NREP increased the S phase, and decreased the G1 phase.Besides that, NREP restrains cell apoptosis in BC cells, which is similar to its role in gastric cancer cell apoptosis [10].These findings help to explain the promotion effect of NREP in BC cell proliferation.
After that, we explored the potential mechanism of NREP in BC.In this work, we confirmed that the expression of NREP could be upregulated by HIF-1α.As we know, HIF-1α is the master regulator of glucose metabolism [38].Additionally, NREP could activate TGF-β1 [22] and Rac1 [25][26][27].TGF-β1 could promote malignant progression of tumor cells by inducing glycolysis [23,24], and inhibition of Rac1 inhibits glycolysis and proliferation of BC cells [28].Therefore, we suspected that NREP, which is upregulated by HIF-1α, may affect glycolysis of BC cells.As we mentioned before, glycolysis is a distinctive cellular metabolic manner in cancer cells exhibiting an increasing rate of glucose uptake and lactic acid fermentation.Therefore, we detected the effect of NREP on glycolysis-related factors.NREP upregulation increased cellular glucose consumption, ATP level, and lactate production.Besides that, NREP promoted the expression of glucose transporters GLUT1, HK2, PKM2, and LDHA in BC cells.Therefore, these results indicated that NREP may promote the growth and metastasis of BC by modulating glycolysis.However, due to the unavailability of the instrument, extracellular acidification rate (ECAR) detection was missed, which is considered a limitation of our current research.
Apart from cell cycle and glycolysis, "DNA replication" pathway was also significantly enriched in KEGG analysis.This result suggested that NREP may also play a carcinogenic role in BC by affecting the factors that related to DNA replication.In addition, based on RNA-seq results, the mRNA expression of the top 10 down-regulated DEGs (PNMT, SMIM10L2B, SNCB, SLC30A3, PLCB2, PTPRN, ADGRE1, GABRD, NTNG2, and BATF3) was validated in NREP up-/down-regulated MDA-MB-468 cells, and the results confirmed that NREP regulates the expression of these genes (Figure S2).Among these genes, SLC30A3 [60], PLCB2 [61], PTPRN [62,63], GABRD [64,65], and BATF3 [66] have been reported were involved in cell proliferation, migration, or metastasis in tumors.PLCB2 is highly expressed in BC and is associated with a poor outcome [67].Thus, these genes may also be involved in NREP mediated malignant phenotypes of BC progression.Of course, more experiments were needed for verification, and these factors and pathways will also be a key direction for our subsequent exploration of NREP.
Taken together, our study showed evidence that NREP could be transcriptionally activated by HIF-1α.In addition, NREP may promote BC tumor growth and metastasis by modulating glycolysis.In summary, the above findings illustrate the possibility that NREP takes great part in BC progression.

Clinical sample collection
BC clinical patients were collected for IHC staining.All specimens were collected according to the provisions of the Declaration of Helsinki, and approved by the Ethics Committee of Zhengzhou University.

Luciferase reporter assay
Dual luciferase reporter gene assay kit (KGAF040, KeyGEN BioTech, China) according to the manufacturer's instruction.The NREP promoter region was inserted into the pGL3-Basic vector (BR014, Fenghui Bio, China) and co-transfected with pRL-TK, siNC, or si-HIF-1α into BC cells.HIF-1α siRNA was synthesized by General Biosystems (China).HIF-1α siRNA: 5ʹ-CAGAAAUGGCCUUGUGAAATT-3ʹ, 5ʹ-UUUCACAAGGCCAUUUCUGTT-3ʹ.After 24 h, BC cells were cultured in normoxia (N) or hypoxia (H) for another 24 h.After transfection, cells were treated with 250 μL lysis buffer.Then, 100 μL firefly luciferase detection reagent and 20 μL sample were added to each well, and the mixture was gently blown before reading.After that, 100 μL renilla luciferase detection reagent was added to each well.Luciferase activity was measured with a multifunctional microplate reader (M200Pro, TECAN, Switzerland).The relative luciferase activity was calculated by firefly/renilla.

Chromatin immunoprecipitation (ChIP)
The ChIP detection kit (WLA106a, Wanlei Bio, China) was used according to the manufacturer's instructions.Cells were treated with 1% formaldehyde to crosslink chromatin associated proteins to DNA.Cell lysates were subjected to ultrasound for 12 sets of 10 s (interval: 30 s) to obtain DNA fragments.Cell lysates were incubated with anti-HIF-1α (20960-1-AP, Proteintech, China), anti-IgG (negative control), and anti-RNA polymerase II (positive control).All the above chromatin supernatants were incubated with 60 μL magnetic protein A/G beads at 4 °C for 1-2 h with rotation.After that, the protein-DNA complexes were reversed and purified by a DNA Gel Extraction Kit (Wanlei Bio, China).DNA products were verified by PCR.Primer for site 2 was forward 5ʹ-TGACTTAGCCCACAGACACT-3ʹ and reverse 5ʹ-AAAGCAGAAGCCATGAATAT-3ʹ, and for site 4 was forward 5ʹ-AACTATC-CACCATCTTG-3ʹ and reverse 5ʹ-ACACCACCTACCACTG-3ʹ.

MTT assay
MTT detection kit (C0009, Beyotime, China) was used in this section.MTT solution (10 μL) was added to each well and the cells were incubated for Y. Ruan et al. 4 h.Formazan (100 μL) was added and continued to incubate until dissolved.The OD value at 570 nm was measured by a microplate reader (800TS, BIOTEK, USA).

Cell cycle
Cell cycle distribution in BC cells was measured by a detection kit (KGA9101, KeyGEN BioTech, China).After centrifugation, cells were fixed in 70% ethanol overnight at 4 °C.Then, cells were washed and resuspended in 500 μL PI/RNase A working solution.Cells were incubated for 30 min in the dark, followed by flow cytometry detection (NovoCyte, Agilent, USA).

Colony formation assay
Cells (400 per dish) were seeded and the colonies could be formed about 15 days later.Cells were fixed with 4% paraformaldehyde (C104188, Aladdin, China) for 25 min at room temperature.Then, Weigert-Giemsa solution (KGA227, KeyGen Biotech, China) was added and stained for 5 min.Colonies containing more than 50 cells were counted under a microscope.Colony formation rate = (number of colonies /number of inoculated cells) ×100%.

EdU staining
EdU imaging detection kit (KGA331, KeyGen Biotech, China) was used in this section.Cells were treated with 10 μM EdU staining solution and incubated for 2 h at 37 °C.Then, cells were fixed with 4% paraformaldehyde for 15 min and 0.5% Triton X-100 for 20 min (0.1 mL) at room temperature.Finally, the click-it solution was added and incubated for 30 min in the dark.After stained with DAPI, images were taken under a fluorescence microscope (IX53, OLYMPUS, Japan).

Tube formation assay in HUVECs
HUVEC cells were cultured in a dedicated medium (PriMed-iCell-002, iCell, China).Matrigel (30 µL; 356234, Corning, USA) was plated in 96-well plates for 2 h at 37 °C in the incubator.Then, HUVECs were seeded, and BC cellderived conditioned medium was added into wells.Then, cells were cultured in an incubator with 5% CO 2 at 37 °C.After 6 h, tubule formation was assessed using microscopy.

Cell apoptosis
Apoptosis in BC cells was measured by Annexin V-FITC/PI double staining apoptosis detection kit (KGA1102, BioTech, China).After centrifugation, cells were resuspended in 500 μL binding buffer.Then, 5 μL AnnexinV-FITC and 5 μL PI were separately added and mixed.Cells were incubated for 15 min in the dark, followed by flow cytometry detection (NovoCyte, Agilent, USA).

Animal experiments
Four weeks old female BALB/c mice (SPF, Cavens, China) were divided into the following groups: Parental, shNC, NREP-sh1, NREP-sh2, Parental, Vector, and NREP, based on the random number table approach.Six mice were included in each group to ensure the accuracy and completeness of the experimental results.The sample size was estimated by power analysis and our experience.No animals were excluded.
For each mouse, 5 × 10 6 cells were subcutaneously injected into the groin.Tumor volumes were measured every 4 days after sizeable tumor formation.At the end of the experiment (24 days later), mice were sacrificed, and tumors were dissected and weighed.For tumor lung metastasis, the mouse was injected with 2.5 × 10 6 cells via the tail vein.The mice were sacrificed after nine weeks, and the lungs were removed and collected for further analysis.

TUNEL staining
Cell apoptosis in mice formed tumor tissues was measured by In Situ Cell Death Detection Kit (12156792910, Roche, Switzerland).Sections (5 μm) were added with 0.1% Triton X-100 (ST795, Beyotime, China) for 8 min at room temperature.Then, the TUNEL working solution was prepared and incubated for 1 h at 37 °C in the dark.DAPI (D106471-5mg, Aladdin, China) was used for counterstaining.Images were taken under a microscope (BX53, OLYMPUS, Japan).

Transwell assay
Cell migration and invasion were measured by transwell assay.For cell migration, the transwell chamber (14341, LABSELECT, China) was placed into a 24-well plate, 800 µL of culture medium containing 10% FBS was added to the lower chamber, and 200 µL of cell suspension was added to the upper chamber, followed by 24 h incubation at 37 °C in the incubator with 5% CO 2 .For cell invasion, transwell chambers were placed in 24-well plates, coated with pre-diluted matrigel gel (356234, Corning, USA), and placed in an incubator at 37 °C for 2 h.After that, a transwell chamber containing matrigel gel was placed into a 24-well plate, 800 µL of culture medium containing 10% FBS was added to the lower chamber, and 200 µL of cell suspension was added to the upper chamber, followed by 48 h incubation.After washing, the chambers were fixed with 4% paraformaldehyde and stained with crystal violet staining solution (0528, Amresco, USA) for 5 min.The cells were counted under a microscope (IX53, OLYMPUS, Japan).Five random fields were selected from each sample, and the average was taken.

RNA-sequencing (RNA-seq)
NREP downregulated MDA-MB-468 cells were cultured and collected for RNA-seq to analyze the downstream genes that may be regulated by NREP.RNA integrity was assessed using the Fragment Analyzer 5400 (Agilent Technologies, CA, USA).DEGs were selected with |Log 2 Foldchange | > 1 and p < 0.05.

Statistical analysis
All evaluators were blinded to the assignments, and all data were used for analysis.Data were shown as mean ± standard deviation (S.D.).Statistical analysis was performed using GraphPad 8.0 statistical software.Normality and homogeneity of variance were tested at first.For comparison between the two groups, an unpaired t-test was employed according to the homoscedasticity of variances, and Mann Whitney test was applied to the non-normally distributed data.For comparisons among multiple groups, one-way ANOVA followed Tukey's multiple comparisons test was applied to the normally distributed and homogeneous data, and Kruskal-Wallis's test was applied to the non-normally distributed data.P values less than 0.05 were considered statistically significant.

Fig. 1
Fig. 1 NREP is upregulated and predicts poor prognosis value in BC patients.A Venn diagram of 22 upregulated differentially expressed genes (DEGs) in nine BC databases from GEO. B Heat map of 22 DEGs in nine BC databases.C Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of 22 DEGs was analyzed by DAVIVD.D Gene ontology (GO) enrichment analysis of 22 DEGs was analyzed by DAVIVD.BP biological process, CC cellular component, MF molecular function.E The mRNA expression level of NREP in nine BC databases.T: tumor tissues; C: normal tissues.F The expression of NREP in the tumor, para-carcinoma, and healthy tissues of BC patients (analyzed by Breast Cancer Gene-Expression Miner v4.9).G Survival curves of BC patients with high or low NREP level (analyzed by Biomedical Informatics Institute).Overall survival (OS), metastasis-free survival (MFS), progression-free interval (PFI), and progression-free survival (PFS).H IHC staining of BC clinical samples with low or high NREP expression.Data were expressed as mean ± SD. *P < 0.05, **P < 0.01.

Fig. 2
Fig. 2 NREP is transcriptional activated by HIF-1α.A The expression of NREP in hypoxia and normoxia BC cells (data from GSE3188 and GSE111259).B The expression of NREP in HIF-1α depleted BC cells (data from GSE3188).C The mRNA expression of NREP in different BC cell lines (SK-BR-3, MCF-7, MDA-MB-468, MDA-MB-231) at normoxia or hypoxia was measured by qPCR.The mRNA expression of HIF-1α (D) and NREP (E) in HIF-1α inhibited MDA-MB-468 or MDA-MB-231 cells were measured by qPCR.F The protein expression and quantification data of NREP and HIF-1α in MDA-MB-468 or MDA-MB-231 cells were measured by western blot.G, H Luciferase reporter assay was used to confirm the binding sites of HIF-1α on NREP promoter.I, J ChIP-PCR assay was used to confirm that HIF-1α could directly bind with the promoter of NREP.Data were expressed as mean ± SD. *P < 0.05, **P < 0.01.

Fig. 3
Fig. 3 NREP promotes BC cell proliferation and angiogenesis in vitro.NREP was down-/up-regulated in MDA-MB-468 and MDA-MB-231 cells.A, B The mRNA expression of NREP in MDA-MB-468 or MDA-MB-231 cells was measured by qPCR.C, D The protein expression of NREP in MDA-MB-468 or MDA-MB-231 cells was measured by western blot.E, F Cell proliferation was analyzed by MTT assay in MDA-MB-468 or MDA-MB-231 cells.G, H The cell cycle was analyzed by flow cytometry in MDA-MB-468 or MDA-MB-231 cells.I, J Colony formation images and quantification data in MDA-MB-468 or MDA-MB-231 cells.K, L EdU staining images (green light) of MDA-MB-468 cells or MDA-MB-231 cells.DAPI: blue.M Cell apoptosis was analyzed by flow cytometry in MDA-MB-468 or MDA-MB-231 cells.N Caspase-3 activity in MDA-MB-468 or MDA-MB-231 cells.O, P Tube formation images and quantification data of HUVECs in tumor-conditioned medium.Data were expressed as mean ± SD. *P < 0.05, **P < 0.01.

Fig. 4
Fig. 4 NREP promotes BC cell migration, invasion, and metastasis.A, B Cell migration was measured by transwell assay in NREP down-/upregulated MDA-MB-231 cells.C, D Cell invasion was measured by transwell assay in NREP down-/up-regulated MDA-MB-231 cells.E The morphological changes of MDA-MB-231 cells were observed under a light microscope.F The protein expression and quantification data of E-cadherin, N-cadherin, SLUG, and MMP9 in NREP down-/up-regulated MDA-MB-231 cells were analyzed by western blot.G, H The mRNA expression of E-cadherin and N-cadherin in NREP down-/up-regulated MDA-MB-231 cells was analyzed by qPCR.I The images of metastasis lesions in mouse lung tissues.J, K H&E staining images and quantification data of metastasis lesions in mouse lung tissues.Data were expressed as mean ± SD. *P < 0.05, **P < 0.01.

Fig. 5
Fig. 5 NREP promotes BC tumor growth in vivo.A Tumor image of xenograft mice injected with NREP-shRNA or -OE infected MDA-MB-468 cells.B Tumor volume changes in the nude mice that were injected with NREP-shRNA or -OE infected MDA-MB-468 cells were analyzed over time.C, D The mRNA expression of NREP and Ki67 in mice formed tumor tissues was analyzed by qPCR.E, F The IHC images and positive cells of NREP and Ki67 in mice formed tumor tissues.G, H TUNEL staining images of mice formed tumor tissues.I IHC images of CD31 and microvessel density of mice formed tumor tissues.J, K The mRNA expression of E-cadherin and N-cadherin in mice formed tumor tissues was analyzed by qPCR.Data were expressed as mean ± SD. *P < 0.05, **P < 0.01.

Fig. 7
Fig. 7 RNA-seq and functional enrichment analyses in NREP silenced BC cells.A, B Principal component analysis (PCA) plot of sh-NC and sh-NREP infected MDA-MB-468 cells.C Hierarchical clustering analysis was performed based on all DEGs.D Volcano map of all DEGs.The red dot stands for upregulated DEGs, blue dot stands for downregulated DEGs.E Selected KEGG enrichment pathway of all DEGs.F The mRNA expression of CCNA2 and CCNE1 in NREP knockdown or overexpressed MDA-MB-468 cells.G The mRNA expression of LDHA, GLUT1, PKM2, and HK2 in NREP suppressed MDA-MB-468 cells (Data from RNA-seq).Data were expressed as mean ± SD. *P < 0.05, **P < 0.01.